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Abstract 

The influence of secondary orientation on 
the elastic response of a [001] -or i ented nickel- 
base single-crystal superalloy, PWA 1480, was 
investigated under mechanical, thermal, and com- 
bined thermal and mechanical loading conditions 
by applying finite element techniques. Elastic 
stress analyses were performed with a commer- 
cially available finite element code. Secondary 
orientation of the single-crystal superalloy was 
offset with respect to the global coordinate 
system in increments from 0° to 90° and stresses 
developed within the single crystal were deter- 
mined for each loading condition. The results 
indicated that the stresses were strongly influ- 
enced by the angular offset between the secondary 
crystal orientation and the global coordinate 
system. The degree of influence was found to 
vary with the type of loading condition (mechani- 
cal, thermal, or combined) imposed on the single- 
crystal superal loy . 

Introduction 

Single crystal (SC) superalloys were iden- 
tified as potential blade materials that war- 
ranted further development for advanced liquid 
propellant rocket engine turbopumps by Chandler 
(1) , After subsequent development, materials 
scientists evaluated the mechanical properties of 
several single crystal superalloys (2,3) to find 
a suitable candidate replacement blade material 
for the turbopumps of the Space Shuttle Main 
Engine (SSME) and advanced rocket engines. A 
well characterized, nickel-base cubic single 
crystal superalloy, PWA 1480, was selected by 
Pratt & Whitney as the blade material for the 
Alternate Turbopump Development (ATD) program for 
the SSME. 

Turbine blades of nickel-base single crys- 
tal superalloys (including PWA 1480) are direc- 
tionally solidified along the low modulus [001] 
crystallographic direction to enhance thermal 


fatigue resistance. The directional solidifica- 
tion process usually generates a secondary crys- 
tallographic direction, [010] that is randomly 
oriented with respect to fixed geometric axes in 
the turbine blade. However, orientation of the 
secondary crystallographic direction can be 
controlled by using a seed crystal during solidi- 
fication (4). Since single crystals exhibit 
anisotropic elastic behavior, the stress-strain 
response and dynamic characteristics of a turbine 
blade that is directionally cast along the [001] 
crystal orientation tend to vary with the orien- 
tation of the secondary crystallographic 
direction within the blade (5-6). 

The main objective of this study was to 
determine the influence of secondary orientation 
on the elastic response of [001] -or i ented, single 
crystal PWA 1480 under mechanical, thermal, and 
combined thermal and mechanical loading condi- 
tions by applying finite element stress analyses. 
A parametric study was conducted in which the 
secondary orientation angle was varied in incre- 
ments of 10° from 0° to 90° for each loading 
condition. The stresses developed within the 
single crystal PWA 1480 superalloy for all the 
cases of the parametric study are presented. 

The turbine blade of a rocket engine turbo- 
pump is schematically represented in Fig. 1. 

This is a hollow SC PWA 1480 turbine blade that 
is directionally cast such that the [001] crys- 
tallographic direction (primary orientation) is 
parallel to the span of the blade. In this par- 
ticular design the hollow core extends below the 
platform into the shank region. In Fig. 1, XYZ 
is the global coordinate system, with Z-axis 
along the span of the blade, Y-axis along the 
chord of the blade, and X-axis along the thick- 
ness of the blade. 

In a cryogenic liquid propellant rocket 
engine turbopump, the blade attachment region may 
be actively cooled. Since the turbine blade air- 
foil is exposed to the hot gas flow path, thermal 
gradients are induced both along the span of the 
blade as well as through the thickness, for the 



full length of the hollow core. These gradients 
may be especially severe through the thickness 
(of the order of 875 °C/cm) , while being milder 
along the span. These gradients are substan- 
tially higher than those in typical aircraft gas 
turbine blades. To examine the role which the 
secondary crystal orientation may play in such 
cases, an analysis of a square plate subjected to 
similar loadings will be employed (Fig. 1). 

A square plate with the dimensions of 25.4 
by 25.4 by 3.2 mm was selected for the parametric 
finite element stress analyses. The crystal 
(X c Y c Z c ) and global (XYZ) coordinate systems for 
the plate are shown in Fig. 1. The global 
coordinate system is obtained by rotating the 
crystal coordinate system about the Z -axis by 
an angle of 6. Since the axes Z and Z c are 
both oriented along the [001] crystallographic 
direction, the angle 6 between the crystal and 
global coordinate systems represents the orien- 
tation of the secondary crystallographic direc- 
tion. The matrix of direction cosines that 
relates the global coordinate system to the 
crystal coordinate system (Fig. 1) is as follows: 


cos # sin # 0 |x 

c 

-sin# cos# 0 Y c (!) 

0 0 1 




x sin 2 # cos 2 # 


D i; = c 12 + 2 < C u - C 12 ~ 2C „> (4) 

x sin 2 # cos 2 # 


ELASTIC STIFFNESS COEFFICIENTS 


Three independent elastic stiffness coeffi- 
cients (C u , C 12 and C 4A ; Eq. (2)) are required to 
describe the elastic behavior of a cubic single 
crystal such as SC PWA 1480 within the crystal- 
lographic coordinate system (7). To describe the 
elastic behavior in the global coordinate system 
the stiffness tensor must be transformed from the 


crystal coordinate system to the global coor- 
dinate system (Fig. 1) according to the law of 
transformation of fourth rank tensors. Lieberman 


and Zirinsky (8) developed a simplified method 
for transforming the stiffness coefficients of 
single crystals between two coordinate systems 
that are related by a matrix of direction cosines 
(Eq. (1)). Their method was used together with 
Eq. (1) to obtain the elastic stiffness coeffi- 
cients of a cubic single crystal in the global 
coordinate system. The stiffness matrix of the 
cubic single crystal within the global coordinate 
system and its relationship to the engineering 
stresses and engineering strains is as follows: 


D is 3 < c n ~ c u -2C 44 ) (sin 3 0cos0 (5) 
- sin#cos 3 #) 



respect i vely. and e . t . e and 7,7, 
r . XX' yy' zz ' , fyz> ' zx* 

7 xy are the engineering normal and shear 
strains, respectively. Thus, Eqs. (2) to (7) 
define the stiffness matrix in a global 
coordinate system (XYZ) that is offset by an 
(arbitrary) angle 0 to the crystal coordinate 
system (X c Y c Z c ; Fig. 1). 
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The three independent elastic stiffness 
coefficients (applicable within the crystal coor- 
dinate system) for SC PWA 1480 at various temper- 
atures are listed in Table 1. These coefficients 
were determined by the Engineering Division of 
Pratt A Whitney, United Technologies Corporation 
under the NASA Lewis Research Center contract 
NAS3-23939. The variations of stiffness coeffi- 
cients D u , D„, D J8 , D and D with the 
secondary orientation angle a, determined trom 
Eqs . (2) to (7) for SC PWA 1480 at 38 °C, are 
shown in Fig. 2. It is evident from this figure 
that the variation in the elastic stiffness 
coefficients has a periodicity of 90° in 9 , and 
therefore, a range of 0 to 90° was selected for 
the secondary orientation angle in the parametric 
study . 

FINITE ELEMENT ANALYSES 

Elastic stress analyses were conducted on 
the plate model shown in Fig. 3 by using the MARC 
finite element structural analysis code (9) . The 
finite element model of the plate contained a 
total of 500 isoparametric eight-noded sol id 
elements with ten elements along the side and 
five elements through the thickness of the plate. 
The boundary conditions used were as follows: 

(1) the node at D was fixed in X, Y, and Z 
directions, (2) nodes along the edge AD were 
fixed in Y and Z directions, (3) nodes along 
the edge DC were fixed in X and Z directions 
and (4) nodes along edges AB and BC were 
fixed in the Z direction. All boundary 
conditions were enforced by specifying the fixed 
boundary option of the MARC code. A total of 
five different loading cases were considered in 
this study. In these cases mechanical, thermal 
and combined mechanical and thermal loads were 
applied separately to the model. The previously 
described boundary conditions were employed in 
all five loading cases. 

The first loading case involved a distrib- 
uted mechanical load of 137.9 MPa along the 
Z-direction (Fig. 4(a)). In the second and third 
loading cases thermal gradients were imposed 
along the Z- and X-direction, respectively 
(Figs. 4(b) and (c)) . Imposition of thermal 
gradients along the Z and X directions was an 
attempt to simulate the thermal gradients that 
occur along the span and through the thickness of 
the actual turbine blade. Finally, the fourth 
and fifth loading cases were obtained by combin- 
ing the first loading case with the second and 
third loading cases, respectively. 

The physical properties of SC PWA 1480 
including mean coefficient of thermal expansion, 
thermal conductivity and specific heat were 
obtained from Ref. 10. For the loading cases 
involving thermal gradients, thermal boundary 
conditions represented by the temperatures in 
Fig. 4 were applied to the finite element model 
considered. These thermal boundary conditions 
may be representative of cryogenic rocket engine 
turbomachinery turbine blade applications. Ther- 
mal boundary conditions were introduced into the 


finite element code enforcing the assigned tem- 
peratures to the specific faces of the plate in 
Fig. 4, while insulating all other faces. Ther- 
mal gradients shown in Fig. 5 were used in the 
second and third loading cases. These data 
represent temperatures at the centroids of the 
elements. The thermal output feature of MARC 
code was used to generate the input subsequent 
stress analysis. 

For each loading case in the parametric 
study, the secondary orientation angle 9 was 
varied from 0° to 90° in increments of 10°. In 
addition, all loading cases were studied for 
6 = 45° to obtain the data at the symmetric loca- 
tion (Fig. 2). The stiffness matrices for the 
various secondary orientations were incorporated 
in the finite element analyses by programming 
Eq. (2) in the MARC user's subroutine HOOKLW (9). 
Variations in the physical and elastic properties 
of the material with temperature were accounted 
for in the analyses. 

The variations in the different components 
of stress with secondary crystal orientation 
angle 6 were documented for three elements shown 
in Fig. 3. These elements were selected to 
capture the extreme as wel I as reference stress 
locations for the five different loading cases 
considered in this study. For example, under the 
second loading case in which a thermal gradient 
was imposed along the Z-di recti on (Fig. 4(b)), 
extreme thermal stresses could be generated in 
faces ABCD and EFGH. Hence, the selection of 
element 295 (Fig. 3) which was near surface EFGH 
for this loading case. Similar reasoning was 
used in selecting element 445 (Fig. 3), for the 
third loading case in which a thermal gradient 
was imposed along the X-direction. Finally, a 
centrally located, element 245 (Fig. 3) was 
selected for the first loading case in which only 
the mechanical loading was imposed. Another con- 
sideration that led to the choice of these spe- 
cific locations was mitigation of the constraint 
effects on the stresses generated within the 
elements. 

For the five loading cases considered in 
this study, the values of the six stress compo- 
nents at the centroid of the selected elements 
are listed in Table 3 for a secondary orientation 
angle of 0 = 0°. In all loading cases the 
stresses developed within the three selected 
elements were well within the elastic regime of 
SC PWA 1480. In loading Case 1, which involved 
only mechanical loading, values of the six stress 
components in the three selected elements were 
very similar. However, in loading Cases 2 and 3 
which involved thermal gradients along the Z- 
and X-axis, respectively, elements 295 (near the 
hot side of a thermal gradient along Z-axis) and 
445 (near the cold side of a thermal gradient 
along X-axis) were subjected to large normal 
compressive and tensile stresses, respectively, 
compared to the normal stresses in element 245 
which was nearly at the center of the mesh. The 
magnitudes of the tensi le normal stresses of 
element 445 in loading Case 3 were larger than 
the magnitudes of compressive normal stresses in 
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element 295 under loading Case 2 because the 
square plate was subjected to a more severe 
thermal gradient along the thickness (X-axis) 
than along the length (Z-axis). Finally, in 
loading Cases 4 and 5 which involved both mech- 
anical and thermal loading, the stresses devel- 
oped within the selected elements were very 
nearly equal to those expected from the super- 
position of stresses developed in the individual 
mechanical and thermal loadings (loading Case 1 
with loading Cases 2 and 3, respectively). These 
results are in agreement with the principle of 
superposition for the linear elastic analysis. 

The influence of secondary orientation on 
the stresses developed within the selected 
elements are shown in Figs. 6 to 8. For each 
element and loading case, the individual stress 
components at a given secondary orientation angle 
6 were normalized with the corresponding stress 
components obtained for 5=0°, i .e. 

(ff W (O,, and (r %y ) e /(T ) and 
so on. "ftie indi^*iaual as well as normalized 
stress components were then plotted against the 
secondary orientation angle 9 for each element 
and loading condition. If there was no influence 
of secondary orientation on the elastic response 
of SC PWA 1480, then all six normalized stress 
components should remain at a fixed value of 
unity as 9 is varied from 0° to 90°. Any 
deviation from unity is an indication that the 
elastic response of PWA 1480 SC is affected by 
the secondary orientation angle 9. The influ- 
ence of secondary orientation on individual and 
normalized stress components under imposed 
mechanical loading (Case 1) is shown in Fig. 6 
for the three selected elements. The variation 
of individual and normalized stress components 
with secondary orientation angle 9 are shown in 
Figs. 7 and 8 for thermal loads (loading Cases 2 
and 3) and combined thermal and mechanical loads 
(loading Cases 4 and 5), respectively. 

The influence of the secondary orientation 
angle on both the individual and normalized 
stress components was dictated by the type (mech- 
anical or thermal or combined) of loading imposed 
on the SC PWA 1480 square plate. For example, no 
significant influence of secondary orientation 
was observed for the individual stress components 
under mechanical loading (Fig. 6). However, in 
cases involving thermal gradients along X- and 
Z-axis the individual normal stress components 
exhibited significantly larger influence of sec- 
ondary orientation than the individual shear 
stresses (Figs. 7 and 8). The normalized stress 
components exhibited a substantial dependence on 
the secondary orientation angle, 9 under all 
five loading cases. In all of the loading cases 
except those involving thermal gradient along the 
X-axis, largest influence of secondary orienta- 
tion was observed for the normalized stress com- 
ponent (j )J{r For the ca ? es * nvo I v i ng 

thermal gradient along the X-axis the normalized 
components (o ) 0 /(ff xx ) (J and (0^/(0, 
exhibited about the same order of variation with 
the secondary orientation angle, 9. 


DISCUSSION 

In this study, the elastic stress compo- 
nents were normalized to determine the influence 
of secondary orientation angle 9 , under differ- 
ent loading conditions. In Figs. 6 to 8 the 
normalized stress components indicate that, on a 
relative basis, all six components of stress are 
influenced by the secondary orientation angle, 9 . 
However, the importance of secondary orientation 
effects cannot be evaluated by consideration of 
the variation in the normalized stress components 
(with 9) alone. The magnitudes of the 
individual components must also be taken into 
consideration for accurate assessment of the 
influence of secondary orientation. For example, 
one normalized stress component, ( T %y ) q/ i r X y) 
exhibited a range as large as ±166 in element 295 
under loading Case 2, in which thermal loading 
was imposed along the Z-axis (Fig. 7 (a) ( i i ) ) . 
However, it is clear from Table 3 that the value 
of r at 6 = 0° is only 0.005 Mpa. For this 
loading condition, multiplying this stress com- 
ponent by a factor of 166 yields 0.830 MPa, a 
relatively small value compared to the other 
normal components of stress within this element 
(Fig. 7(a) (i)). By comparison, the normal ized 
stress component (VyJg/ ( a yy )^ increased only 
from 1.0 to 1.4 as 7 v was increased from 0° 
to 90° under thermal loading along the X-axis 
(Fig. 7 (b) ( i i ) ) • However, this represented an 
increase in o from 204 2 MPa at 0 = 0° to 
286.3 MPa at ¥ = 45° (Fig. 7(b) (ii)). Thus, the 
variations in the normalized stress components 
with 9, while indicative of the influence of 
secondary orientation on the elastic response of 
SC PWA 1480, alone are not sufficient to assess 
the importance of secondary orientation effects. 
The variation of actual magnitudes of the indi- 
vidual components of stress also must be consid- 
ered to determine the quantitative significance 
of the secondary orientation effects. 

Among the mechanical, thermal and combined 
loading cases considered in this study, the case 
of a thermal gradient along the X-axis (through 
the thickness of the square plate) produced large 
stresses (Table 3; In Case 3, element 445 is 
under biaxial tension). This is due to the fact 
that the thermal gradient through the thickness 
of the plate is more severe than the thermal 
gradient along the length of the square plate 
(Fig. 5). In loading case 5, the mechanical load 
which was superimposed on the thermal load of 
Case 3, nearly doubled the o zz component of 
stress within the square plate (Table 3; Case 5, 
e I ement 445) . 

The state of maximum biaxial tension 
occurred for loading Case 5 at a secondary 
orientation angle 9 = 45° (Table 3, element 445 
and Fig. 8(b) (i)). This is to be expected 
because in this instance the stiffness coeffi- 
cient D.j (which reaches a maximum at 9 = 45°, 
Fig. 2) is coupled with a steep thermal gradient 
through the thickness of the square plate and a 
mechanical load along the Z-axis. Thus, large 
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thermal gradients perpendicular to the primary 
orientation [001] of the turbine blade may be 
expected to produce very high stresses, if they 
occur near the secondary orientation of 0 - 45°. 
However, these stresses can be minimized by 
controlling the secondary orientation so that 
thermal gradients perpendicular to [001] occur 
near 0=0°. 

Some normalized stress components shown in 
Figs. 6 to 8 exhibit a lack of symmetry about the 
secondary orientation angle, 6 = 45°. For exam- 
ple, the normalized stress component $/ ( a xx ) 0 

is not symmetric about 6 = 45°. Similar 
behavior is also exhibited by the normalized 
stress components, ) and (r ) g/ (r ) 0 

i" Fig. 6(c) (ii) and & ' 

The main reason for lack or symmetry about 
0 = 45° is the coupling between <r„„ and a., 
and 7 x as well as between r 
e in X the single crystal stiffness 
(Eq. (2)). The coupling elastic coefficients 
D a and D 2a change sign about 6 = 45® (Fig. 2) 
which creates a lack of symmetry in the stress 
components. A small amount of observed lack of 
symmetry may also have been caused by numerical 
round off errors. This may especially be the 
reason for the observed asymmetry in 
(T )^/(r ) , because r has no coupling terms 
in the single crystal stiffness matrix. 

In this study the influence of secondary 
orientation was studied under purely elastic 
conditions and steady state thermal gradients. 

In actual service, turbine blades are subjected 
to severe thermal transients due to start up and 
shutdown, which can produce localized plasticity 
and creep in the turbine blades. The influence 
of secondary orientation on the stresses produced 
within the actual blade under conditions of plas- 
ticity and creep will be different from the 
influence determined under fully elastic condi- 
tions. To determine the influence of secondary 
orientation on the stresses developed within the 
turbine blades, nonlinear, transient thermal 
structural analyses must be conducted with appro- 
priate boundary conditions. 


CONCLUSIONS 


A parametric study was conducted with 
elastic finite element analyses to determine the 
influence of secondary orientation angle on the 
elastic behavior of a [001] -or iented nickel-base 
single crystal superalloy, PWA 1480 under 
mechanical, thermal and combined mechanical and 
thermal loads. The secondary orientation angle 
was varied from 0° to 90° and the stresses 
developed within a square plate of SC PWA 1480 
under different loading conditions were computed. 
The following conclusions were drawn from the 
parametric study. 

1. Under the loading conditions studied, 
the normalized stress components showed a sub- 
stantial variation with the secondary orientation 
angle. However, the normalized stress components 
should be considered along with the magnitudes of 
the individual stress components to determine the 


significance of secondary orientation in 
engineering design. 

2, The type of loading (mechanical, ther- 
mal, or combined) imposed on the square plate 
dictated the variation in the individual stress 
components with the secondary orientation angle. 

3. For the same temperature boundary condi- 
tions, thermal gradients through the thickness of 
the square plate produced higher stresses than 
thermal gradients along the length (primary 
orientation of [001]) of the square plate. Thi~ 
was due to a higher thermal gradient through t.<-j 
thickness of the plate. If it is required to 
impose thermal gradients perpend i cu I ar to the 
primary [001] orientation, a secondary 
orientation of 0° tends to minimize the thermal 
stresses developed within the plate, whereas a 
secondary orientation of 45° tends to produce 
much larger thermal stresses. 
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TABLE 1.- INDEPENDENT ELASTIC 
STIFFNESS COEFFICIENTS 


FOR SC PWA 1480 


Temperature, 

°C 

c n> 

GPa 

C !2< 

GPa 

^ 44 ' 

GPa 

-18 

252 

163 

131 

38 

250 

163 

129 

93 

248 

161 

128 

149 

246 

160 

126 

204 

244 

159 

124 

260 

242 

158 

123 

316 

240 

157 

121 

371 

238 

157 

119 

427 

235 

156 

117 

482 

233 

154 

115 

538 

230 

154 

113 

593 

228 

152 

111 

649 

225 

152 

109 

704 

223 

151 

107 

760 

219 

150 

105 

816 

217 

150 

102 

871 

213 

149 

100 

927 

210 

148 

97 

982 

206 

146 

95 

1038 

201 

145 

92 

1093 

197 

143 

88 

1149 

192 

143 

85 

1204 

186 

142 

81 


TABLE 2. - PHYSICAL PROPERTIES OF SC PWA 1480 


[From Ref. (9).] 


Temperature, 

Therma 1 

— — 

Therma 1 

Spec i f i c 

°C 

coeff icient 

conduct i v i ty , 

heat, 


of expansion, 

k, 

^ i 


<*, 

W/m °C 

kJ/kg °C 


m/m/°C 



93 

0.0000107 

9.2 

0.406 

204 

.0000111 

10.8 

.440 

316 

.0000115 

12.4 

.465 

427 

.0000119 

13.7 

.477 

538 

.0000124 

15.3 

.481 

649 

.0000129 

16.9 

.490 

760 

.0000134 

18.5 

.511 

871 

.0000140 

20.0 

.569 

982 

.0000147 

21.5 

.687 

1093 

.0000157 

22.8 

.917 


TABLE 3. - STRESSES DEVELOPED IN SINGLE CRYSTAL PWA 1480 PLATE UNDER MECHANICAL, 


THERMAL AND COMBINED LOADING CASES SECONDARY ORIENTATION ANGLE, 6 = 0° 


sms 

Type of 

Element 


°yy' 


T xy> 

V 


II 

loading 


MPa 

MPa 

MPa 

MPa 

MPa 

MPa 

i 

Mechanical 

245 

0.655 

1.310 

136.9 

0.001 

0.316 

—6.047 


(Z-axis) 

295 

.012 

.053 

137.9 

.001 

-.003 

-.006 


445 

-.090 

1.234 

136.6 

.005 

.385 

-.182 

2 

therma I 

245 

.043 

-15.78 

19.98 

-0.017 

-5.475 

-.296 


(Z-axis) 

295 

-19.44 

-20.34 

-1.88 

.005 

.043 

-0.034 

3 

Therma 1 

245 

.034 

3.020 

6.847 

1.496 

-.270 

-12.41 


(X-axis) 

445 

.993 

204.2 

130.0 

1.214 

15.56 

-3.744 

4 

Meehan ica 1 

245 

.697 

-14.47 

156.9 

.018 

-5.157 

-.343 


and thermal 
(Z-axia) 

295 

-19.43 

—20.29 

136.0 

.006 

.040 

-.041 


Mechanics 1 

245 

.687 

4.330 

143.7 

1.496 

.046 

12.46 

5 

(Z-axis) 
and thermal 

445 

.903 

205.4 

266.4 

1.220 

15.95 

-3.923 


(X-axis) 




- 
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Z, Z c [001] 
H 



Z [001 ] 




XYZ : Global co-ordinate system. 
X C Y C Z C : Crystal co-ordinate system. 


Figure 1. — Schematics of SSME turbine blade and the square plate considered for the analysis. 







Temperature ("C) 




(c) Thermal gradient along X-axis. 

Figure 4. -Loading cases imposed on the square plate. 



Distance along Z-axis (mm) Distance along X-axis (mm) 

(a) Along the height of the square plate (Z-axis). (b) Through the thickness of the square plate (X-axis). 
Figure 5. -Temperature distributions calculated from the heat transfer analyses. 
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(i) Individual stress components. 

(a) Element 245. 


(ii) Normalized stress components. 



(i) Individual stress components. (ii) Normalized stress components. 

(b) Element 295. 



ANGLE OF SECONDARY ORIENTATION. 0 (Ifc-prees) ANGLE OF SECONDARY ORIENTATION. 0 (Degrees) 


(i) Individual stress components. (ii) Normalized stress components. 

(c) Element 445. 

Figure 6. -Influence of secondary orientation 9 on individual and normalized stress components under imposed 
mechanical loading. 
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(i) Individual stress components. 
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(a) Thermal gradient along Z-axis; Element 295. 


(ii) Normalized stress components. 
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(i) Individual stress components. (ii) Normalized stress components, 

(b) Thermal gradient along X-axis; Element 445. 

Figure 7.-Variation of normalized stress components with secondary orientation angle 9 for Ihermal loads. 




(i) Individual stress components. (ii) Normalized stress components, 

(a) Thermal gradient along Z-axis; Element 295. 




(i) Individual stress components, (ii) Normalized stress components. 

(b)Thermal gradient along X-axis; Element 445. 

Figure 8.-Variation of normalized stress components with secondary orientation angle 9 for combined mechanical 
and thermal loads. 
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